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Amendments to the Claims: 

This listing of claims will replace all prior versions, and listings, of claims in the 

application: 

1. (Presently Amended) A method for reconstructing an image of a scattering medium, 
comprising: 

directing energy into the scattering medium at a source location on the scattering 
medium; 

measuring the energy emerging from the scattering medium at a detector location on the 
scattering medium; 

selecting an initial guess of internal properties of the scattering medium; 

predicting the energy emerging from the scattering medium using an i .ntegro-di f ferential 
equation of radiati ve transfe r that contains at least a streaming term, and, an integral term that 
accounts for photons being scattered from all directions into a specific direction, and wherein the 
prediction is a function of the initial guess; 

generating an objective function that is proportional to the difference between bas e d on a 
comparison of the prediction given by the integro-differential equation of radiative transfer and 
the actual with the measurement; 

generating a gradient of the objective function by a method of adjoint differentiation; 

modifying the initial guess of the properties of the scattering medium based on the 
gradient of the objective function; and 

generating an image representation of the internal properties of the scattering medium. 
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2. (Previously Presented) The method according to claim 1, further comprising repeating 
the predicting of the energy emerging from the scattering medium based on the modified initial 
guess, generating the objective function and modifying the initial guess, until at least one of a 
predetermined number of repetitions has occurred or the objective function reaches a 
predetermined threshold. 

3. (Previously Presented) The method according to claim 1, wherein the prediction depends 
on boundary conditions. 

4. (Presently Amended) The method according to claim 3, wherein the boundary conditions 
account for a refractive index mismatch at an interface between the medium and at least one 
detectors or source. 

5. (Original) The method according to claim 1, wherein the prediction comprises an 
iterative process producing intermediate results. 

6. (Original) The method according to claim 5, wherein the intermediate results of the 
prediction are stored. 

7. (Original) The method according to claim 6, wherein generating the gradient of the 
objective function by adjoint differences uses the intermediate results of the prediction. 

8. (Original) The method according to claim 7, wherein generating the gradient comprises 
stepping backward through the intermediate results of the prediction. 

9. (Presently Amended) The method according to claim 1, wherein the integro-differential 
equation of radiati ve transfer is time independent. 
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10. (Presently Amended) The method according to claim 9, wherein the time independent 
integro-differential equation of radiative transfer is: 



where 4*(r,<2>) is the radiance at the spatial position r directed into a unit solid angle a> 9 S(r,w) is 
the energy directed into the medium at spatial position r into a unit solid angle a> 9 ju s is the 
scattering coefficient, ju a is the absorption coefficient and p{co,co t ) is the scattering phase 
function. 

11. (Original) The method according to claim 10, wherein the scattering phase function is: 



where 0is the angle between the two unit solid angles co and a> 
and g is the anisotropy factor. 

12. (Presently Amended) The method according to claim 1, wherein the integro-differential 
equation of radiative transfer is time dependent. 

13. (Presently Amended) The method according to claim 12, wherein the time dependent 
integro-differential equation of radiative transfer is: 



p(cos0) = 



2(l + g 2 -2gcos^) 



3/2 



1 fflF(r,co,Q 
c dt 




where ^P(r,^y,t) is the radiance at the spatial position r directed into 



a unit solid angle co, S(r,w, t) is the energy directed into the 



medium at spatial position r into a unit solid angle a> 9 ju s is the 
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scattering coefficient, /j a is the absorption coefficient and p{co,co r ) 
is the scattering phase function. 

14. (Original ) The method according to claim 13, wherein the scattering phase function is: 

P(oos0) = ~~ 7 

2(1 + g- ~2gcos60 

where 0 is the angle between the two unit solid angles co and co \ 
and g is the anisotropy factor. 

15. (Previously Presented) The method according to claim 1, wherein the properties include 
at least one of a scattering coefficient, an absorption coefficient, an anisotropy factor, or a 
scattering phase function. 

16. (Previously Presented) The method according to claim 1, wherein the objective function 
is a comparison of the predicted energy and the measured energy. 

17. (Previously Presented) The method according to claim 1, wherein the objective function 
is based on the sum of the differences between the predicted energy and the measured energy for 
each source detector pair, wherein a source detector pair is formed between each source location 
and each detector location. 

18. (Original) The method according to claim 1, wherein the objective function is: 

where Mi represents the actual measurements and the Pi represents the predicted 
measurements for each source defector pair i, m is the number of source detector pairs, where a 
source detector pairs is formed between each source location and each detector location. 
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19. (Original) The method according to claim 1, further comprising minimizing the objective 
function, 

20. (Original) The method according to claim 19, wherein minimizing the objective function 
includes a one dimensional line search. 

2 1 . (Original) The method according to claim 20, wherein the one dimensional line search is 
performed along a direction of the gradient of the objective function. 

22. (Original) The method according to claim 20, wherein the one dimensional line search is 
performed along a gradient-dependent direction. 

23. (Original) The method according to claim 1, wherein the energy comprises near infra-red 
energy. 

24. (Original) The method according to claim 1, wherein the scattering medium contains 
regions wherein the scattering coefficients are not substantially greater than the absorption 
coefficients. 

25. (Original) The method according to claim 1, wherein the scattering medium contains a 
low scattering region embedded in a high scattering region. 

26. (Original) The method according to claim 1, wherein the predicted energy is determined 
using finite element methods. 

27. (Original) The method according to claim 1, wherein the predicted energy is determined 
using finite difference methods. 
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28. (Presently Amended) A method for imaging the spatial optical properties of tissue, 
comprising: 

(a) directing energy into the scattering medium at a source location on the tissue; 

(b) measuring the energy emerging from the scattering medium at a detector location 
on the tissue; 

(c) selecting an initial guess of the spatial optical properties of the tissue; 

(d) predicting the energy emerging from the tissue using an integro-differential 
equation of radiative transfer in an iterative process, wherein the prediction is a function of the 
initial guess and a refraction index mismatch at a boundary of the tissue, and the iterative process 
generates a plurality of intermediate predictions; 

(e) generating an objective function based on a comparison of the prediction with the 
measured energy emerging from the scattering medium; 

(f) generating a gradient of the objective function by adjoint differentiation; 

(g) modifying the initial guess of the spatial properties of the tissue based on the 
gradient of the objective function; 

(h) repeating steps (d) through (g) until at least one of a threshold of modifications to 
the initial guess is reached or the objective function reaches a threshold; and 

(j) generating an image representation of the spatial optical properties of the tissue. 

29. (Presently Amended) A system for reconstructing an image of a scattering medium, 
comprising: 

a source for directing energy into the scattering medium at a source location on the 
scattering medium; 
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a detector for measuring the energy emerging from the scattering medium at a detector 
location on the scattering medium; 

an initial guess of internal properties of the scattering medium; 

means for predicting the energy emerging from the scattering medium using an integro- 
differential equation of radiative transfer, wherein the prediction is a function of the initial guess; 

means for generating an objective function based on a comparison of the prediction with 
the measurement; 

means for generating a gradient of the objective function by a method of adjoint 
differentiation; 

means for modifying the initial guess of the properties of the scattering medium based on 
the gradient of the objective function; and 

means for generating an image representation of the internal properties of the scattering 
medium. 

30. (Previously Presented) The system according to claim 29, further comprising means for 
repeating the predicting of the energy emerging from the scattering medium based on the 
modified initial guess, generating the objective function and modifying the initial guess, until at 
least one of a predetermined number of repetitions has occurred or the objective function reaches 
a predetermined threshold. 

31. (Previously Presented) The system according to claim 29, wherein the prediction 
depends on boundary conditions. 

32. (Previously Presented) The system according to claim 31, wherein the boundary 
conditions account for a refractive mismatch at an interface between the medium and at least one 
detector or source. 
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33. (Previously Presented) The system according to claim 29, wherein the prediction 
comprises an iterative process producing intermediate results. 

34. (Original) The system according to claim 33, wherein the intermediate results of the 
prediction are stored. 

35. (Previously Presented) The system according to claim 34, wherein generating the 
gradient of the objective function by adjoint differentiation uses the intermediate results of the 
prediction. 

36. (Original) The system according to claim 35, wherein generating the gradient comprises 
stepping backward through the intermediate results of the prediction. 

37. (Presently Amended) The system according to claim 29, wherein the integro-differential 
equation of radiative transfer is time independent. 

38. (Presently Amended) The system according to claim 37, wherein the time independent 
integro-differential equation of radiative transfer is: 

am?(r 9 a>) + (//„ + fi s )W(r 9 a>) = W S(r,a>) + ju s ^ p(co,oJ)^{r^)dcd 

where ^(r,co) is the radiance at the spatial position r directed into a unit solid angle co, S(r,w) is 
the energy directed into the medium at spatial position r into a unit solid angle co, /u s is the 
scattering coefficient, ju a is the absorption coefficient and p{(o,af) is the scattering phase 
function. 

39. (Original) The system according to claim 38, wherein the scattering phase function is; 
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p(cos6>) = ■ 



1-g 



2 



2(l + g 2 -2gcos6>) 3/2 
where 0is the angle between the two unit solid angles co and co and g is the 
anisotropy factor. 

40. (Presently Amended) The system according to claim 29, wherein the mtegro-differential 
equation of radiative transfer is time dependent. 

4L (Presently Amended) The system according to claim 40, wherein the time dependent 
integro-differential equation of radiative transfer is: 



1 d¥(r,<D,0 



'■S(r,aj)-®-Sr¥(r,aj)--(ju a + ju s )*¥(r 9 a>,t) + ju s p(co,a) x y¥(r,a)\t)dco } 



c dt 

where ^(r^t) is the radiance at the spatial position r directed into a unit solid angle co, S(r,w, t) 
is the energy directed into the medium at spatial position r into a unit solid angle a), ju s is the 
scattering coefficient, ju a is the absorption coefficient and pico.co 1 ) is the scattering phase 
function. 

42. (Original) The system according to claim 41, wherein the scattering phase function is: 



l-# 2 

p(co$6) = 



v3/2 



20 + g 2 ~-2gcos#) 3 
where #is the angle between the two unit solid angles a> and co \ and g is the 
anisotropy factor. 

43. (Previously Presented) The system according to claim 29, wherein the properties include 
at least one of a scattering coefficient, an absorption coefficient, an anisotropy factor, or a 
scattering phase function. 
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44. (Previously Presented) The system according to claim 29, wherein the objective function 
is a comparison of the predicted energy and the measured energy. 

45. (Previously Presented) The system according to claim 29, wherein the objective function 
is based on the sum of the differences between the predicted energy and the measured energy for 
each source detector pair, wherein a source detector pair is formed between each source location 
and each detector location. 

46. (Previously Presented) The system according to claim 29, wherein the objective function 
is; 

<p = -l(P i -M i ) 2 

where Mj represents the actual measurements and Pi represents the predicted 
measurements for each source detector pair, m is the number of source detector pairs, where a 
source detector pairs is formed between each source location and each detector location. 

47. (Previously Presented) The system according to claim 29, further comprising minimizing 
the objective function. 

48. (Original) The system according to claim 47, wherein minimizing the objective function 
includes a one dimensional line search. 

49. (Original) The system according to claim 48, wherein the one dimensional line search is 
performed along a direction of the gradient of the objective function, 

50. (Original) The system according to claim 49, wherein the one dimensional line search is 
performed along a gradient-dependent direction. 

1025284 vl 



Serial No. 09/767,230 -19- Docket No. 0887-4150US1 

51. (Original) The system according to claim 50, wherein the energy comprises near infra- 
red energy. 

52. (Previously Presented) The system according to claim 29, wherein the scattering medium 
contains regions wherein the scattering coefficients are not substantially greater than the 
absorption coefficients. 

53. (Previously Presented) The system according to claim 29, wherein the scattering medium 
contains a low scattering region embedded in a high scattering region. 

54. (Previously Presented) The system according to claim 29, wherein the predicted energy 
is determined using finite element methods. 

55. (Previously Presented) The system according to claim 29, wherein the predicted energy 
is determined using finite difference methods. 

56. (Presently Amended) A system for imaging the spatial distribution of optical properties 
of tissue, comprising: 

(a) a source for directing energy into the scattering medium at a source location on 
the tissue; 

(b) a detector for measuring the energy emerging from the scattering medium at a 
detector location on the tissue; 

(c) an initial guess of spatial optical properties of the tissue; 

(d) means for predicting the energy emerging from the tissue using an integro- 
differential equation of radiati ve transfer in an iterative process, wherein the prediction is a 
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function of the initial guess and a refraction index mismatch at a boundary of the tissue, and the 
iterative process generates a plurality of intermediate predictions; 

(e) means for generating an objective function based on a comparison of the 
prediction with the measured energy emerging from the scattering medium; 

(f) means for generating a gradient of the objective function by adjoint 
differentiation; 

(g) means for modifying the initial guess of the spatial properties of the tissue based 
on the gradient of the objective function; 

(h) means for repeating steps (d) through (g) until at least one of a threshold of 
modifications to the initial guess is reached or the objective function reaches a threshold; and 

(j) means for generating an image representation of the spatial optical properties of 
the tissue. 

57. (Presently Amended) Computer executable software code stored on a computer readable 
medium, the code for reconstructing an image of a scattering medium, comprising: 

code to direct energy into the scattering medium at a source location on the scattering 
medium; 

code to measure the energy emerging from the scattering medium at a detector location 
on the scattering medium; 

code to receive an initial guess of internal properties of the scattering medium; 

code to predict the energy emerging from the scattering medium using an integro- 
differential equation of radiative transfer, wherein the prediction is a function of the initial guess; 

code to generate an objective function based on a comparison of the prediction with the 
measurement; 
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code to generate a gradient of the objective function by a method of adjoint 
differentiation; 

code to modify the initial guess of the properties of the scattering medium based on the 
gradient of the objective function; and 

code to generate an image representation of the internal properties of the scattering 
medium. 

58. (Presently Amended) Computer executable software code stored on a computer readable 
medium, the code for imaging the spatial distribution of optical properties of tissue, comprising: 

(a) code to direct energy into the scattering medium at a source location on the tissue; 

(b) code to measure the energy emerging from the scattering medium at a detector 
location on the tissue; 

(c) code to receive an initial guess of spatial optical properties of the tissue; 

(d) code to predict the energy emerging from the tissue using an inte gro - di fferenti al 
equation of radiative transfer in an iterative process, wherein the prediction is a function of the 
initial guess and a refraction index mismatch at a boundary of the tissue, and the iterative process 
generates a plurality of intermediate predictions; 

(e) code to generate an objective function based on a comparison of the prediction with 
the measured energy emerging from the scattering medium; 

(f) code to generate a gradient of the objective function by adjoint differentiation; 

(g) code to modify the initial guess of the spatial properties of the tissue based on the 
gradient of the objective function; 

(h) code to repeat steps (d) through (g) until at least one of a threshold of modifications to 
the initial guess is reached or the objective function reaches a threshold; and 
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(j) code to generate an image representation of the spatial optical properties of the tissue. 

59. (Presently Amended) A computer readable medium having computer executable 
software code stored thereon, the code for reconstructing an image of a scattering medium, 
comprising: 

code to direct energy into the scattering medium at a source location on the scattering 
medium; 

code to measure the energy emerging from the scattering medium at a detector location 
on the scattering medium; 

code to receive an initial guess of internal properties of the scattering medium; 

code to predict the energy emerging from the scattering medium using an integro- 
differential equation of radiati ve transfer, wherein the prediction is a function of the initial guess; 

code to generate an objective function based on a comparison of the prediction with the 
measurement; 

code to generate a gradient of the objective function by a method of adjoint 
differentiation; 

code to modify the initial guess of the properties of the scattering medium based on the 
gradient of the objective function; and 

code to generate an image representation of the internal properties of the scattering 
medium. 

60. (Presently Amended) A computer readable medium having computer executable 
software code stored thereon, the code for imaging the spatial distribution of optical properties of 
tissue, comprising: 

(a) code to direct energy into the scattering medium at a source location on the tissue; 
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(b) code to measure the energy emerging from the scattering medium at a detector 
location on the tissue; 

(c) code to receive an initial guess of spatial optical properties of the tissue; 

(d) code to predict the energy emerging from the tissue using an integro-differential 
equation of radiative transfer in an iterative process, wherein the prediction is a function of the 
initial guess and a refraction index mismatch at a boundary of the tissue, and the iterative process 
generates a plurality of intermediate predictions; 

(e) code to generate an objective function based on a comparison of the prediction with 
the measured energy emerging from the scattering medium; 

(f) code to generate a gradient of the objective function by adjoint differentiation; 

(g) code to modify the initial guess of the spatial properties of the tissue based on the 
gradient of the objective function; 

(h) code to repeat steps (d) through (g) until at least one of a threshold of modifications to 
the initial guess is reached or the objective function reaches a threshold; and 

(j) code to generate an image representation of the spatial optical properties of the tissue. 



1025284 vl 



